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The Marginalized (5-GLMB Filter 

Claudio Fantacci, Ba-Tuong Vo, Francesco Papi and Ba-Ngu Vo 


Abstract 

The multi-target Bayes filter proposed by Mahler is a principled solution to recursive Bayesian tracking based on 
RFS or FISST. The 5-GLMB filter is an exact closed form solution to the multi-target Bayes recursion which yields 
joint state and label or trajectory estimates in the presence of clutter, missed detections and association uncertainty. 

Due to presence of explicit data associations in the 5-GLMB filter, the number of components in the posterior grows 
without bound in time. In this work we propose an efficient approximation to the 5-GLMB filter which preserves both 
the PHD and cardinality distribution of the labeled posterior. This approximation also facilitates efficient multi-sensor 
tracking with detection-based measurements. Simulation results are presented to verify the proposed approach. 

Index Terms 

RFS, FISST, 5-GLMB filter, LMB filter, PHD 

I. Introduction 

Multi-target filtering/tracking involves the simultaneous estimation of the number of targets along with their 
states, based on a sequence of noisy measurements such as radar or sonar waveforms |[T]. To reduce complexity and 
facilitate tractability, the sensor waveforms are typically processed into a sequence of detections. The key challenges 
in multi-target filtering/tracking thus include detection uncertainty, clutter, and data association uncertainty. To date, 
three major approaches to multi-target tracking/filtering have emerged as the main solution paradigms. These are, 
Multiple Hypotheses Tracking (MHT), 0-0> Joint Probabilistic Data Association (JPDA) Q, Q, and Random 
Finite Set (RFS) Q. 

The RFS or Finite Set STatistics (FISST) approach pioneered by Mahler provides principled recursive Bayesian 
formulation of the multi-target filtering/tracking problem. The essence of the RFS approach is the modeling of the 
collection of target states and measurements, referred to as the multi-target state and multi-target measurement, as 
finite set valued random variables (T), 0. The centerpiece of the RFS approach is the Bayes multi-target filter (T), 
which recursively propagates the filtering density of the multi-target state forward in time. The PHD Q, Q, CPHD 
03 and cardinality-balanced and labeled Multi-Bernoulli filters pT| , fl^ are tractable approximations to the 
Bayes multi-target filter which are synonymous with the RFS framework. Their tractability however largely hinges 
on the approximate form for the posterior which cannot accommodate statistical dependencies between targets. 
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The Bayes multi-target filter is also a (multi-target) tracker when target identities or labels are incorporated into 
individual target states. In HD, |T4), the notion of labeled RFSs is introduced to address target trajectories and their 
uniqueness. The key results include conjugate priors that are closed under the Chapman-Kolmogorov equation, and 
an analytic solution to the Bayes multi-target tracking filter known as the ^-Generalized Labeled Multi-Bernoulli 


(J-GLMB) filter 1151. With detection based measurements, the computational complexity in the i5-GLMB filter is 
mainly due to the presence of explicit data associations. For certain applications such as tracking with multiple 
sensors, partially observable measurements or decentralized estimation, the application of a 5-GLMB filter may 
not be possible due to limited computational resources. Thus cheaper approximations to the (5-GLMB filter are of 
practical significance in multi-target tracking. 

In this paper we present a new approximation to the 5-GLMB filter. Our result is based on the approximation 
proposed in fib) where it was shown that the GLMB distribution can be used to construct a principled approximation 
to an arbitrary labeled RFS density that matches the PHD and the cardinality distribution. We refer to the resultant 
filter as a Marginalized 5-GLMB (Mi5-GLMB) filter since it can be interpreted as a marginalization over the 
data associations. The proposed filter is consequently computationally cheaper than the (5-GLMB filter while 
still preserving key summary statistics of the multi-target posterior. Importantly the M5-GLMB filter facilitates 
tractable multi-sensor multi-target tracking. Unlike PHD/CPHD and Multi-Bernoulli based filters, the proposed 
approximation accommodates statistical dependence between targets. We also present an alternative derivation of 


the LMB filter proposed in |12| based on the newly proposed M5-GLMB filter. Simulations results verify the 
proposed approximation. 


II. Background and problem formulation 

This section briefly presents background material on multi-object filtering and labeled RFS, which form the basis 
for the formulation of our multi-target tracking problem. 


A. Multi-object Estimation 

Suppose that at time k, there are N{k) object states Xk.i, ■ ■ ■ each taking values in a state space 

X. In the random finite set (RFS) framework, the multi-object state at time k is represented by the finite set 
Xk = , Xk.Nik)}, and the multi-object state space is the space of all finite subsets of X, denoted as -F(T’). 

An RFS is simply a random variable that take values the space TiX') that does not inherit the usual Euclidean notion 
of integration and density. Mahler’s Finite Set Statistics (FISST) provides powerful yet practical mathematical tools 
for dealing with RFSs ||T] based on a notion of integration/density that is consistent with point process theory m- 

Let TTk{‘\Zk) denote the multi-target posterior density at time k, and iTk+iiki’) denote the multi-target prediction 
density to time k H- 1 (formally Tr^j-) and 'Xk+i\k{ ) should be written respectively as 7rfc(-|Zo,..., Zk-i, Zk), and 
7rfc_|_i|fe(-|Zo) • ■ •! Zk), but for simplicity the dependence on past measurements is omitted). Then, the multi-target 
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Bayes recursion propagates Trk{ ) in time 0, 0, according to the following update and prediction 

s gk{Zk\Xk)Tri~^i~_i{Xk) 

'^k\-^k\^k) — ~~p 5 

J gk{Zk\X)TTk\k-iiX)SX 

7Tk+i\kiXk) = j fk\k-i{Xk\X)Trk-iiX)SX, (2) 

where fk\k-i is the multi-object transition density to time fc + 1, is the multi-object likelihood function at time 
k, and the integral is a set integral defined for any function f : X{X) ^ 'R hy 

j f{X)dX = J f{{xi,...,Xi})d{xi,...,Xi). (3) 

2 = 0 

An analytic solution to the multi-object Bayes filter for labeled states and track estimation from the multi-object 
filtering density was given in GZl- 

B. Labeled RFS 

To perform tracking in the RFS framework we use the label RFS model that incorporates a unique label in the 
object’s state vector to identify its trajectory 0. In this model, the single-object state space A” is a Cartesian product 
XxL, where X is the feature/kinematic space and L is the (discrete) label space. A finite subset set X of XxL 
has distinct labels if and only if X and its labels {£ : {x, i) G X} have the same cardinality. An RFS on XxL with 
distinct labels is called a labeled RFS di- 

For the rest of the paper, we use the standard inner product notation (/, p) = / f{x)g{x)dx, and multi-object 
exponential notation — 11 ^(®)’ where /i is a real-valued function, with /i® = 1 by convention. We denote a 

generalization of the Kroneker delta and the inclusion function that take arbitrary arguments such as sets, vectors, 
etc, by 

1, if X = Y 


6y{X) 4 

li-(X) 4 


0, otherwise 


1 , if X c y 

0, otherwise 

We also write lyjx) in place of lY({a;}) when X = {x}. Single-object states are represented by lowercase letters, 
e.g. X, X while multi-object states are represented by uppercase letters, e.g. X, X, symbols for labeled states and 
their distributions are bolded to distinguish them from unlabeled ones, e.g. x, X, tt, etc, spaces are represented by 
blackboard bold e.g. X, Z, L, etc. 

An important class of labeled RFS is the generalized labeled multi-Bernoulli (GLMB) family p7) , which is the 
basis of an analytic solution to the Bayes multi-object filter fTS] . Under the standard multi-object measurement 
model, the GLMB is a conjugate prior that is also closed under the Chapman-Kolmogorov equation. If we start with 
a GLMB initial prior, then the multi-object prediction and posterior densities at any time are also GLMB densities. 
Let £ : XxL —>■ L be the projection £((x,f)) = £, and A(X) A(5|^| (|£(X)|) denote the distinct label indicator. 
A GLMB is a labeled RFS on XxL distributed according to 


7r(X) = A(X)^u;(^)(£(X)) 


cGC 




n X 


( 4 ) 
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where C is a discrete index set, (L) and satisfy: 


= 1, (5) 

LCL cGC 

J p^‘^\x,t)dx = 1 . ( 6 ) 

The GLMB density Q can be interpreted as a mixture of multi-object exponentials. Each term in Q consists of a 
weight w^‘^\£{'X.)) that depends only on the labels of X, and a multi-object exponential that depends on 

the entire X. The PHD (or intensity function) of the unlabeled version of generalized labeled multi-Bernoulli RES 
is given by 

cGC t&u LCL 

The Labeled Multi-Bernoulli (LMB) family is a special case of the GLMB family with one term: 


p^^\x,t) = p^^\x) 

= n(i--■">) 

.^gm leL 

where ,p^^'>)}£^M, M C L, is a given set of parameters with representing the existence probability of 

track i, and p^^'> the probability density of the kinematic state of track £ given its existence | [l^ . Note that the index 
space C has only one element, in which case the (c) superscript is not needed. The LMB family is the basis of the 
LMB filter, an effective approximation of the Bayes multi-target tracking filter, which is highly parallelizable and 
capable of tracking large number of targets GD- The LMB filter, however, is an approximation of the Bayes multi¬ 
target tracking filter which only preserves the unlabeled PHD of the multi-target posterior m- The information 
lost (e.g. the cardinality distribution of multi-target posterior and the approximate construction of the individual 
tracks) can lead to poor performance in cases of low observability and/or signal-to-noise ratio (SNR) which will 
be demonstrated in Section m 


C. The 6-GLMB Filter 


An efficient approach to multi-target tracking was presented in | [T7| using a special form of the GLMB distribution 
in eq. Q called (5-GLMB, i.e. 


7r(X) = A(X) Y, 

(/,5)GJG(L)xH 



(7) 


= A(X) ^ 5z(>C(X))^u;(^'«) 



( 8 ) 


{GH 

The i5-GLMB density naturally arises in multi-target tracking problems when using the standard detection based 
measurement model. In the following we briefly recall the prediction and update steps for the J-GLMB filter, 
additional details can be found in HD- To ensure distinct labels we assign each target an ordered pair of integers 
£ = {k, i), where k is the time of birth and i is a unique index to distinguish targets born at the same time. The label 
space for targets born at time fc-|-1 is denoted as L^+i, and a target born at time k + 1, has state x G XxL/^+i. The 
label space for targets up to time k + 1 (i.e. including those born prior to fc -f 1), denoted as Lo:fe+i, is constructed 
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recursively by Lo:fc+i = Lo:fe UL^+i (note that Lo:fc and L^+i are disjoint). A multi-object state X at time fc + 1, 
is a finite subset of A = XxLo:/c+i. 

Suppose that at time k, there are Nk objects with states Xfc i,..., each taking values in the (labeled) 

state space X x Lo:fc, and measurements i, ..., s^^h taking values in an observation space Z. The 

multi-object state and multi-object observation, at time k, Q, 0 are, respectively, the finite sets 


II 


(9) 

Yfc = {j/fc.i,. 

■ • : yk,Mk } ■ 

(10) 


The 5-GLMB filter recursively propagates a 5-GLMB posterior density forward in time according to the following 
Bayesian update and prediction 


'TT’/c \^k') — 

9k {^k A:—1 (^fc) 

(11) 

J 5fc(Xfc|X)7rfe|fc_i(X)5x’ 


J /fe|fe_i(Xfc|X)7rfc_i(X)<5X, 

(12) 

which is the labeled counterpart of the Bayesian 

recursion ([T])-Q. 



1) S-GLMB Prediction: Given the current multi-object state X', each state G X' either continues to 

exist at the next time step with probability Ps{x',P) and evolves to a new state {x,£) with probability density 
fk+i\kix\x',f)6i{£'), or dies with probability 1 — Ps{x' ,£'). Note that the label of the objects is preserved in the 
transition, only the kinematic part of state changes. Assuming that X has distinct labels and that conditional on 
X, the transition of the kinematic states are mutually independent, then the set W of surviving objects at the next 
time is a labeled multi-Bernoulli RFS dz) 

fs(W|X) = A(W) A(X) l£(x) (£(W)) [$ (W; , (13) 

where 

<i>(W;x',/)= E 5i.{£)Ps{x',£')f {x\x', £')+[! {£')]{! -Ps{x',£')) . (14) 

(x,^)GW 

The A(X) in ( [T3| l ensures that only X with distinct labels are considered. The set of new objects bom at the next 
time step is distributed according to 

fB(Y) = A(Y)u;b(/:(Y)) [pb]^ (15) 


The birth density fg is defined on X x L^+i and fB(Y) = 0 if Y contains any element y with £(y) ^ L 


fe-i-i- 


The birth model (15i covers both labeled Poisson and labeled multi-Bernoulli. The multi-object state at the next 
time X is the superposition of surviving objects and new born objects, i.e. X = W U Y. Since the label spaces L 
and B are disjoint, the labeled birth objects and surviving objects are independent. Thus the multi-target transition 
density turns out to be the product of the transition density (fO]) and the density of new objects ([T5]l 


f (X|X') = fs(X n (X X L) |X') fs (X - (X X L)) . 

Additional details can be found in GD- 


(16) 
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If the current multi-object prior density is a (5-GLMB of the form then the multi-object prediction density is 
a (5-GLMB given by 


7rfe+i|fc(X) = A(X) 


E 




(/,^)GJ^(Lo;fc+l)xE 


'M) 

i^k+l\k 


n X 


where 


(17) 


^k+l\k 

= UlB(/-Lo:fc)w5^(/nLo:fc), 

(18) 

Pf+l\k«) 

= <0:kWPs'' O + <k+iWPB{x,£) , 

(19) 

P^i\x',£) 

(Psi;e)fk+i\k{x\;e),pi^H;e)) 


VsH^) 

(20) 

4\£) 

= {Ps{;e),pi^\;£)) , 

(21) 



(22) 


J Q]Lo:fc 



2) 6-GLMB Update: The standard multi-object observation model is described as follows. For a given multi-object 
state X, each state x € X is either detected with probability Pjj (x) and generates a point z with likelihood g{z\x.), 
or missed with probability 1 — Pd (x), i.e. x generates a Bernoulli RFS with parameter {Pd{:x-), p(-|x)). Assuming 
that conditional on X these Bernoulli RFSs are independent, then the set W C Z of detected points (non-clutter 
measurements) is a multi-Bernoulli RFS with parameter set {(P£i(x), p(-|x)): x G X}. The set F C Z of false 
observations (or clutter), assumed independent of the detected points, is modeled by a Poisson RFS with intensity 
function k{-). The multi-object observation Z is the superposition of the detected points and false observations, i.e. 
Z = W U Y, and the multi-target likelihood can be derived as shown in HI)- Assuming that, conditional on X, 
detections are independent, and that clutter is independent of the detections, the multi-object likelihood is given by 


y] , (23) 

6/ee(z;(x)) 

where ©(/) is the set of mappings 6* : / —>■ {0,1,..., M}, such that 9{i) = 9{i') > 0 implies i = i' , and 

^j:Y{x,(.\9) = I <99(1)) . (24) 

[ l-PD(,x,e), if 0 (£)=O 

Note that an association map 9 specifies which tracks generated which measurements, i.e. track ^ generates 
measurement ye(t) € Y, with undetected tracks assigned to 0. The condition “9{i) = 9{i') > 0 implies i = T”, 
means that a track can generate at most one measurement, and a measurement can be assigned to at most one track, 
at one time instant. 

If the current multi-object prediction density is a (5-GLMB of the form Q, then the multi-object posterior density 
is a i5-GLMB given by 


TTfc (X|Z) 


A(x) y] y] 




n X 




(25) 
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where 0(/) denotes the subset of the current maps with domain I, and 

n I 


0^ ^fe+l|fe 


(26) 


P<^ (- 1 ^) =- im77. -■ 

Vz W 

Notice that the new association maps 9 can be added (stacked) to their respective association histories ^ in order 
to have again the more compact form (0 for the updated i5-GLMB (25 i. 


III. The Marginalized J-GLMB Filter 

In this section we present a new solution for recursive multi-target tracking based on the GLMB approximation 
technique presented in flh) . The resultant filter is called the Marginalized (5-GLMB (M5-GLMB) filter since the 
result can be interpreted as performing a marginalization with respect to the association histories. We present two 
important justifications for the new algorithm, namely, a computationally efficient approximation of the Bayes 
optimal ^-GLMB filter which directly facilitates multi-sensor updates, and a theoretical result showing that the 
proposed approximation matches exactly the (labeled) PHD and the cardinality distribution of the filtering density. 
Furthermore we show a connection with the LMB filter by presenting an alternative derivation of the LMB filter 
based on extracting individual tracks from the M(5-GLMB filter. 


A. Marginalized 5-GLMB Approximation 

One of the main factors contributing to the computational complexity of the (5-GLMB filter fTS] is the exponential 
growth of the number of hypotheses in the update of the prior ( [26| which gives rise to the an explicit sum over 
an association history variable. Moreover, in multi-sensor scenarios the number of association histories is further 
increased due to successive update steps (as detailed in subsection |III-B2 1 . The idea behind the proposed M(5- 
GLMB filter is to construct a principled GLMB approximation ff (•) to the posterior density tt (•) which results in 
a marginalization over the association histories thereby drastically reducing the number of components required to 
represent the posterior or filtering density. 


Definition 1. A Marginalized (5-GLMB density ff corresponding to the (5-GLMB density tt in (j^ is a probability 
density of the form 


fr(X) = A(X) ^ ^,(/:(X))u;(^) 

/GJr-(L) 



(27) 


where 


wG) — ^ , 

f63 

p^P\x,t} = . 

^ 5GH 


(28) 

(29) 
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Proposition 1. The Marginalized S-GLMB density tt in (27i-(29i preserves both PHD and cardinality distribution 
of the original 6-GLMB density tt in ([^. 

Proof: We apply the result in Proposition 2 of m which can be used to calculate the parameters of the 
marginalized (5-GLMB density. Notice that the result in m applies to any labeled RFS density and our first 
step is to rewrite the ^-GLMB density 0 in the general form for a labeled RFS density specified in p6| , i.e. 
7r(X) = w(>C(X))p(X) where 

.. ,£n}) = j 7r({(a;i,fi),..., ... ,a;„) 


/GJ='(L) 




ies /gjg(l) 


and 




^l)) ■ • ■ ) (^ri) £n)}) 

W{{ii,. . . ,in}) 




w 


iPi) 


JA) 






/GJ='(L) 






(30) 

Applying Proposition 2 of flh) , the parameters and p^^'^ for the M(5-GLMB approximation that match the 
cardinality and PHD are 

(L) = 6i{L)w{I) = 6i{L) E 


and 


p''G(^x,i) = l/(f)p/-{£}(a;,f) 

= l/(^) J Pi{ix,£),{xi,ii),... ,{xj,£j)})d{xi,... ,Xj) 


(31) 


where we enumerate / — {£} = {£i,... ,£j}. Substituting the expression (30i in (31 1 we have 
p(-f) (x, £) = li{£)A{{{x, £), (xi, £i), ..., {xj,£j)})- 


1 


w({£,£i,...,£j}) 


E Sj{{£,£u---,e,})J2 


w 


(J.?) 


JGJ='(L) 






{{x,e),{xi,ei),...,Gj,ej)} 


d{xi,.. .,Xj) 


= l/(^)A({(a;,£), {xiji ),..., (xj,£j)}) -— E E u>^'^'^'>p^^'>{x,£) 


and noting that I = {£, £i,..., £j} it follows that only one term in the sum over J is non-zero thus giving 

(x,£) = l/(f)A({(a;, £),{xu£i),..., {xj,£j)}) ^ E (a ;,£) 


E' 
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Consequently, the M5-GLMB approximation is given by 


#(X)= ^ 

/GJ='(L) 




n X 


= A(X) ^ <57(/:(X))^u; 

/GJG(L) jgh 


( 1,0 


X 


eeJ 


= A(X) ^ <57(/:(X))r 


M) 


.(0 


X 


/GJ='(L) 


where 






(a;,£) = ^ 


E 




= l/(^)E 




B. M6-GLMB Recursion 

The M(5-GLMB density can be exploited to construct an efficient recursive multi-object tracking filter by cal¬ 
culating the M(5-GLMB approximation step after the ^-GLMB update, and predicting forward in time using the 
i5-GLMB prediction. 


1) MS-GLMB Prediction: Given an updated density of the form (27i the Mi5-GLMB prediction step turns out 
to be 


where 


7r,+i|,(X) = A(X) ^ {L (X)) 


n X 


/GJ^(Lo;fc + l) 


(32) 


‘^k+l\k 

= wb{I- l-o-.k)w^PiI Gho,k), 

( 33 ) 

Pkhlki^^^) 

= ^TLo:kWPs\x',£) + lTL^^,{e)pB{x,£) , 

( 34 ) 

P^p{x',i) 

{Psi;£)fk+i\kix\;£),pi'\;£)) 

Vs\^) 

( 35 ) 

Vs\^) 

= {Ps{;£),pi'\;£)) , 

( 36 ) 

Ws\l) 


( 37 ) 





which is exactly the (5-GLMB prediction step ( [T7| with no association histories from previous time step, i.e. ^ = 0, 
and with the convention of having the superscript (/) instead of (^) due to the marginalization (28l-(29l. 
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Remark 1. The number of components computed after the M5-GLMB prediction step (32 1 is 

IJ^(Lo:fe+i)|. On the other hand, the number of components after the (5-GLMB prediction (17l 

is |J-(Lo:fc+i) X S| for and |S| forp^^|^|^. Notice that the number of weights of the M(5-GLMB is 

substantially lower than the of (5-GLMB. As for the number of location PDFs it is worth noticing 

that the growth rate of the association histories ^ G S is super-exponential with time 115|, 0. while the growth 
rate of the cardinality of J^(Lo:fe+i) is by far more restrained. 


The use of the M(5-GLMB approximation further reduces the number of hypotheses in the posterior density while 
preserving the PHD and cardinality distribution 0- Moreover, the M(5-GLMB is in a form that it is suitable for 
efficient and tractable information fusion (i.e. multi-sensor processing) which will be shown in the next subsection. 

2 ) Multi-Sensor MS-GLMB Update: Consider now a multi-sensor setting in which the sensors (indexed with s) 
convey all the measurement sets to a central fusion node. Assuming that such a measurement sets taken by the 
sensors are conditionally independent on the states, the multi-object Bayesian filtering update 0 can be naturally 
extended as follows: 

n5fc(^l|Xfe)7rfc|fc_i(Xfc) 

7r(Xfc) = -. (38) 

/ Y[gUZl\X)7Tkik-imSX 


where is the multi-object likelihood of sensor s. Thus, at each time instant k, the M(5-GLMB update step (391 


(and equivalently for the (5-GLMB update step (25 i) is sequentially repeated exploiting the measurement sets Z^ 
provided by the sensors. 

Let us now focus on the single update step to be carried out for each sensor s. If the current multi-object prior 
density is a M(5-GLMB of the form ([27]l, then the multi-object posterior density is a (5-GLMB given by 


TTfe (X|Z) = A(X) ^ ^ wi^’'^\z)SiiC{X)) 

/Gjr-(L) eee{i) 


Pk 


{■\z) 


(39) 


where ©(/) (see (23l) denotes the subset of the current maps with domain / , and 


Ufi) U) 

^k ot 




p''k''H-\z) = 


^zix,e; 9) 
Pd{x,£) g{z0(^e)\x,£) 


l-PDix,£), 



(40) 

) 

(41) 


(42) 

if 0{£) > 0 

if 9{£) = 0 

(43) 


Using now (28 i-(29i, the M5-GLMB density corresponding to the (5-GLMB density in (39 1 is a probability density 
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of the form ( [27) t with 


,U) 


= E ^ 

660 ( 7 ) 






(44) 


(45) 


660(7) 


The M(5-GLMB density provided by (44i-(45i preserves both PHD and cardinality distribution of the original S- 
GLMB density. Summing up, at each time instant k, (39 i-(|45]l have to be carried sequentially for each sensor s to 
evaluate 

Remark 2. Each hypothesis I & F (L) generates a set of |0(f)| new measurement-to-track association maps for 
the (5-GLMB posterior. The number of components stored/computed after the M5-GLMB update 

step (39 1 is |-E(L)| • X] 7 gj='(l) I®(^)I' On the other hand, the number of hypotheses after the 


i5-GLMB update (25 1 is |J^(L) x .n] • |0(^)| for 

conclusions along the lines of Remark [T] hold. 


and 


• E76.f(l) |0(^)I for Pi+l|fc- 


The same 


Remark 3. After the marginalization procedure (44i-(45i only |Jr(L)| hypotheses are retained, as all the new 
contributions provided by the association maps |0(/)| are aggregated in a single component. Notice that |Jr(L)| is 
the exact same number of hypotheses produced during the prediction step (|32|) (see Remark [^. Thus, the prediction 


step (39 1 sets the upper bound of the total hypotheses that will be retained after each full M<5-GLMB step. 


From Remark and the MJ-GLMB is preferable over the (5-GLMB in terms of stored information and 
computational burden, since the number of remaining hypotheses after each sensor update step in ( [38] l is always 
set to |J^(L)|. Note that this does not apply to the J-GLMB due to the super-exponential growth as reported in 
Remark [T] This is an important property of the M(5-GLMB since it yields a principled approximation which greatly 
decreases the need of pruning hypotheses w.r.t. the i5-GLMB flS] . In fact, pruning in the 5-GLMB might lead to poor 
performance in multi-sensor scenarios with low SNR (e.g. high clutter intensity, low probability of detection, etc.) 
and limited storage/computational capabilities. For instance, this may happen if a subset of the sensors do not detect 
one or more targets and hypotheses associated to the true tracks are removed due to pruning. Furthermore, from a 
mathematical viewpoint, pruning between corrections generally produces a less informative and order-independent 


approximation to the posterior distribution in eq. (38l. 


C. M6-GLMB Filter Implementation 

The i5-GLMB hlter implementation [TS] applies directly to M(5-GLMB. For a linear Gaussian multi-target model 
it is assumed that l) the single target transition density, likelihood and birth intensity are assumed to be Gaussian; 
It) survival and detection probabilities are constants; III) each single target density is represented as a Gaussian 
mixture. The corresponding Gaussian mixture predicted and updated densities are computed using the standard 


Gaussian mixture update and prediction formulas based on the Kalman filter |15|. In the case of having non¬ 
linear single target transition density and/or likelihood, one can resort to the well known Extended or Unscented 
Kalman Filters p4), |251. On the other hand, for non-linear non-Gaussian multi-target models (with state dependent 
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survival and detection probabilities), each single target density can be represented by a set of weighted particles. 
The corresponding predicted and updated densities are computed by the standard particle (or Sequential Monte 
Carlo) filter 


D. Connection with the LMB Filter 

The LMB filter introduced in | [T2) is a single component approximation to a (5-GLMB density that matches the 
unlabeled PHD. In this subsection we show an alternative derivation of the LMB approximation first proposed in 
through a connection with the M(5-GLMB approximation. Recall that a LMB density is uniquely parameterized 
by a set of existence probabilities and corresponding track densities 

= A{X)w{C{X))p^ (46) 

where 

=n (i(47) 

zgl 

p{x,t) = p^^\x) (48) 


In the following we show that by extracting individual tracks from the M(5-GLMB approximation, we can obtain 
the same expressions for the existence probabilities and state densities originally proposed for the LMB filter in 

& 


Pr*(£(X) 9f) = ;^u;(L) 

L3l 

(49) 

^ li{£)w{I) 

/GJ7(L) 

(50) 

/GJ7(L) 

(51) 

(/,?)eJF(L)xH 

(52) 

Pr*(X9(x,f)) 

Pr*(/:(X) 9 £) 

(53) 

J 7r{{x,e)CX)5X 

j.{i) 

(54) 


where the notation for the numerator in ( |53| l is defined as per |[T] eq. 11.111], while the numerator of ( |54) i follows 
from |[T] eq. 11.112]. Notice that the numerator is precisely the PHD v corresponding to tt, which by Proposition 
2 of p6) exactly matches the PHD v corresponding to tt. Using the results in | |T6| , it can be verified that 

v{x,£) = ^ li{£)w^^^p^^\x,i) (55) 

/GJ7(L) 
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and consequently 




1 

7*(^) 


/G^(L) 

li{£)w^^’^'*p^^\x,£) 

(/.«)GJ=-(L)xH 


(56) 

(57) 


Notice that, however, the property of matching the labeled PHD of the 5-GLMB does not hold for the LMB 


filter, as shown in 112 Section III], due to the imposed multi-Bernoulli structure for the cardinality distribution. 


IV. Numerical Results 

To assess performance of the proposed Marginalized i5-GLMB (M5-GLMB), a 2-dimensional multi-object tracking 
scenario is considered over a surveillance area of 50 x 50 [km^]. Two sensor sets are used to represent scenarios 
with different observability capabilities. In particular: l) a single Radar in the middle of the surveillance region is 
used as it guarantee observability; ll) a set of 3 range-only (Time Of Arrival, TOA), deployed as shown in Fig. 
[T] are used as they do not guarantee observability individually, but information from different sensors need to be 
combined to achieve it. 

The scenario consists of 5 targets as depicted in Fig. For the sake of comparison, the MJ-GLMB is also 


xio'^ xIO'^ 



Fig. 1. Network with 3 TOA sensors. Fig. 2. Target trajectories considered in the simulation experiment. 

The stait/end point for each trajectory is denoted, respectively, by 
•\B. The indicates a rendezvous point. 

compared with the J-GLMB (J-GLMB) fTS) , p7) and LMB (LMB) fT^ filters. The three tracking filters are 
implemented using Gaussian Mixtures to represent their predicted and updated densities 0, |Tg. Due to the non 
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linearity of the sensors, the Unscented Kalman Filter (UKF) p2| is exploited to update means and covariances of 
the Gaussian components. 

The kinematic object state is denoted by x = Px, Py, Py]^, i-C- the planar position and velocity. The motion 

of objects is modeled according to the Nearly-Constant Velocity (NCV) model |Tg-0: 


Xt+l = 


where = 5 and the sampling interval is Tg = 5 [s]. 

The Radar has the following measurement function: 

AiPx - + j {Py - y^)] 


1 

T 

X s 

0 

0 



1 

1 

2^s 

0 

0 

0 

1 

0 

0 

Xt+Wt, 

Q = A 

1 7~'3 
2 ^ s 

s 

0 

0 

0 

0 

1 

Ts 



0 

0 

1 7^4 
4^s 

1 7~'3 

2 s 

0 

0 

0 

1 



0 

0 

1 7^3 
2 -^ s 

T2 

s 


h{x) = 


(58) 


t/ {Px - X^f + {py - 

where (x’’, y^) represents the known position of the Radar and its measurement noise is a^adar = [1 [°]: 100 
The measurement functions of the 3 TOA of Fig. are: 

Hx) = \J{px- xA + {py - y^f , (59) 

where (x^,y^) represents the known position of sensor (indexed with) s. The standard deviation of the TOA 
measurement noise is taken as aroA = 100 [m]. 

The clutter is characterized by a Poisson process with parameter Ac = 15. The probability of target detection is 
Pd = 0.85. 

In the considered scenario, targets pass through the surveillance area with partial prior information for target birth 
locations. Accordingly, a 10-component LMB RFS ttb = 1*^^ been hypothesized for the birth 

process. Table gives detailed summary of such components. Due to the partial prior information on the object 


TABLE I 

Components of the LMB RFS birth process at a given time k. 


rW = 0.09 
p'g’ix) =N'[x-, 

Pb = diag(l0®, 10"*, 10®, 10"*) 


Label 

(k, 1) 

(k, 2) 

(fc. 3) 


[0, 0, 40000, 0]^ 

[0, 0, 25000, 0]^ 

[0, 0, 5000, 0[^ 

Label 

(k, 4) 

(k, 5) 

(k, 6) 

W 

[5000, 0, 0, 0]^ 

[25000, 0, 0, 0)^ 

[36000, 0, 0, 0[^ 


Label 

(k, 7) 

(k, 8) 

(0 

[50000, 0, 15000, 0)^ 

[50000, 0, 40000, 0)^ 

Label 

(k, 9) 

(fc, 10) 

W 

[40000, 0, 50000, 0)^ 

[10000, 0, 50000, 0)^ 
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birth locations, some of the LMB components cover a state space region where there is no birth. Therefore, clutter 
measurements are more prone to generate false targets. 


Multi-target tracking performance is evaluated in terms of the Optimal SubPattern Analysis (OSPA) metric 1231 
with Euclidean distance, p = 2, and cutoff c = 600 [to]. The reported metric is averaged over 100 Monte Carlo trials 
for the same target trajectories but different, independently generated, clutter and measurement noise realizations. 
The duration of each simulation trial is fixed to 1000 [s] (200 samples). 

The three tracking filters are coupled with the parallel CPHD look ahead strategy described in | |T5| , p7) . The 
CPHD |Tg filter. 


A. Scenario 1: Radar 

Figs. W and display the statistics (mean and standard deviation) of the estimated number of targets obtained, 
respectively, with the MJ-GLMB, the 5-GLMB and the LMB. As it can be seen, all the algorithms estimate the 
target cardinality accurately, with no substantial differences. This result indicates that, in the presence of a single 
sensor guaranteeing observability, the approximations made by both MJ-GLMB and LMB are not critical in that 
they provide performance comparable to the 5-GLMB with the advantage of a cheaper computational burden and 
reduced storage requirements. Note that the problems introduced by the rendezvous point (e.g. merged or lost tracks) 
are correctly tackled by all the algorithms. 

Fig. shows the OSPA distance of the algorithms. Note again that, in agreement with the estimated cardinality 
distributions, the OSPA distances are nearly identical. 


B. Scenario 2: 3 TOA 

Figs. Izli and display the statistics (mean and standard deviation) of the estimated number of targets obtained, 
respectively, with the M5-GLMB, the J-GLMB and the LMB. The Mi5-GLMB and the (5-GLMB tracking filters 
estimate the target cardinality accurately, while the LMB exhibits poor performance and higher standard deviation 
due to losing some tracks when 4 or 5 targets are jointly present in the surveillance area. It is worth noticing that 
the M(5-GLMB performs as nearly as identical to the J-GLMB and that the problems introduced by the rendezvous 
point are again correctly tackled. 


Fig. 10 shows the OSPA distance. Note that the OSPA of the Mi5-GLMB is close to the one of (5-GLMB, while 
the LMB shows an overall higher error in agreement with the cardinality error due to losing tracks. 
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Fig. 4. Cardinality statistics for (5-GLMB tracking filter using 1 Radar. 
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Fig. 6. OSPA distance (c = 600 [m], p = 2) using 1 Radar. 
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Fig. 8 
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Fig. 10. OSPA distance (c = 600 [m], p = 2) using 3 TOA. 
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V. Conclusion and future work 

This paper has proposed a novel approximation to the J-GLMB filter with standard point detection measurements. 
The result is based on a principled GLMB approximation to the labeled RFS posterior that matches exactly 
the posterior PHD and cardinality distribution. The proposed approximation can be interpreted as performing 
a marginalization with respect to the association histories arising from the J-GLMB filter. The key advantage 
of the new filter lies in the reduced growth rate of the number of new components generated at each filtering 
step. In particular, the approximation (or marginalization) step performed after each update is guaranteed to 
reduce the number of generated components which normally arise from multiple measurement-to-track association 
maps. Typically, the proposed M(5-GLMB filter requires much less computation and storage especially in multi¬ 
sensor scenarios compared to the J-GLMB filter. Furthermore the proposed M5-GLMB filter inherits the same 
implementation strategies and parallelizability of the J-GLMB filter. A connection and alternative derivation of the 
LMB filter is also provided. Future works will consider distributed estimation with the MJ-GLMB filter. 
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